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In this work we apply two different real-space renormalization-group (RSRG) approaches to the 
anisotropic antiferromagnetic spin-1/2 Heisenberg model on the square lattice. Our calculations 
allow for an approximate evaluation of the T vs. A phase-diagram: the results suggest the existence 
of a critical value of A > 0, at which the critical temperature goes to zero, and the presence of 
reentrant behavior on the critical line between the ordered and disordered phases. This whole 
critical line is found to belong to the same universality class as the Ising model. Our results are 
in accordance with previous RSRG approaches but not with numerical simulations and spin-wave 
calculations. 
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I. INTRODUCTION 

Two-dimensional antiferromagnctism plays an impor- 
tant role in the description of La2Cu04-based high- 
temperature superconductors In these, the spin fluc- 
tuations in the CuC>2 planes are well described by the 
spin-1/2 antiferromagnetic Heisenberg model on a square 
lattice. This, together with the theoretical discussion 
made by Anderson Q], has raised the interest on this 
model, which has been a challenge for decades. 

The Hamiltonian of the anisotropic Heisenberg model 
reads: 

- [3H = K ^2 [(1 - A) (SfSf + SVS?) + SfSfi , (1) 

where K = J/kT is the dimensionless exchange param- 
eter, with K < 0(> 0) for the antiferromagnetic (ferro- 
magnetic) model, k is the Boltzmann constant, T is the 
temperature, the sum is over pairs of nearest-neighbor 
spins i and j, Sf is the A*^ 1 component of the spin-1/2 
Pauli operator on site I, A is the anisotropy parameter, 
and (3 = 1/kT. Note that A = describes the isotropic 
Heisenberg model and for A = 1 we regain the Ising 
model. 

For J > (ferromagnetic interactions), the ground 
state is well known and the Mermin- Wagner theorem ex- 
cludes long-range order at finite temperatures in two di- 
mensions for A = ||. On the other hand, for < A < 1 
there is an easy-axis and the symmetry is the same as for 



the Ising model; therefore, long-range order is possible. 
Indeed, the universality class for all models with values 
of A is the same as for the Ising model. A schematic 
temperature vs. anisotropy phase-diagram for the fer- 
romagnetic model is presented in Fig|l], where all the 
features discussed above are exhibited. 
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FIG. 1. Schematic plot of the temperature vs. anisotropy 
phase-diagram for the two-dimensional anisotropic ferromag- 
netic Heisenberg model, where O stands for the ferromagnetic 
(ordered) phase and D stands for the paramagnetic (disor- 
dered) phase. The arrow indicates that the transitions for 
A / belong to the same universality class as for the Ising 
model (A = 1). Note that the critical temperature line goes 
to zero in the limit of the isotropic Heisenberg model (A = 0). 

For many classical (and some quantum) systems on 
bipartite lattices, there is a mapping of the ferromag- 
netic model onto the antiferromagnetic one. Neverthe- 
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less, there is no such mapping for the Heisenberg model; 
therefore, the ground state of the anisotropic antiferro- 
magnetic Heisenberg (AAH) model is not obtained from 
its ferromagnetic couterpart by flipping all spins in a 
given sub-lattice. In fact, the ground state for J < is 
not even exactly known and has been a matter of debate 
for a long time. However, long-range order was proven 
to be present at T = for A > A c = 0.40 j|, whereas 
the same authors, using extra assumptions, calculated a 
lower value of A c , namely, A c = 0.09 [0. 

Previous RSRG procedures have been able to calcu- 
late the approximate phase-diagram for the AAH on the 
square lattice [^|J^]. The former reference uses a hierar- 
chical lattice to approximate the square one, and per- 
forms a partial trace over internal degrees of freedom, in 
a manner introduced in Ref. . This approach was the 
extension of the Niemeijer-van Leeuwen method to quan- 
tum spin systems. On the other hand, Ref. applies 
the so-called mean-field renormalization group ideas. Al- 
though the approximate scaling transformations calcu- 
lated in Refs. [^) and || are different, the results obtained 
are qualitatively the same, and we will comment on them 
when our results are discussed. Nevertheless, one should 
bear in mind that RSRG procedures on quantum systems 
do not have the same firm basis as on classical mod- 
els. The reason is the non-commutivity aspects of the 
Hamiltonian and, in as concerns Ref. [§, the necessity 
of introducing symmetry-breaking fields which are cho- 
sen according to the ground state of the classical Ising 
model. 

Results from spin-wave theory || on the isotropic 
Heisenberg model give the same value for the critical 
temperature of the ferromagnetic and antiferromagnetic 
systems ]9j , which is in direct contradiction to rigorous 
results |jlO| and to some approximate calculations HH]. 

Therefore, the form of the phase-diagram of the AAH 
model on the square lattice is far from a settled question 
and more work is desirable to study its critical properties. 
We thus apply two different RSRG procedures to evalu- 
ate the approximate phase-diagram of this system. The 
ferromagnetic Heisenberg case was also studied, in order 
to compare our results with previous ones for a model for 
which the critical behavior is well known. 

The remainder of this paper is organized as follows. In 
section II we outline the formalism we used, in section III 
we present results, and in the last section we summarize 
our main conclusions. 



II. FORMALISM 

In this section we outline the finite-size scaling RG 
(FSSRG) procedure, since it is somewhat new in the lit- 
erature and, to the best of our knowledge, has never been 
used in the study of quantum systems. This method was 
proposed some time ago [O and has been successfully 



applied to classical systems (both static and dynamic 
properties have been studied), like Ising, Potts, or Blume- 
Capel models The second is just the generalization 
of the usual bond-moving Migdal-Kadanoff |l3| approxi- 
mation to antiferromagnetic quantum systems. 

The finite-size scaling assumption is that, near the crit- 
ical region, thermodynamic quantities have the following 
form 111: 



(2) 



where b is some arbitrary scaling factor, L is the lin- 
ear dimension of the lattice, T is a scaling function, 
e = \T — T c \ (T c is the critical temperature), and v is 
the critical exponent of the correlation length, such that 
£ ~ e~ v for T close to T c and in the thermodynamic 
limit. The exponent ip is the anomalous dimension of 
the thermodynamic quantity for the magnetization 
M, tp — —[3/v, while for the magnetic susceptibility \i 
if) — 7/iA 

The idea behind the FSSRG is to construct quantities 
which have a zero anomalous dimension, ijj = 0, such 
that: 



Qu{e') = Q v {b 1 ' v e) = Q L {e). 



(3) 



As seen, these quantities will have the same value at 
T = T c , no matter what the lattice sizes L or L' are 
(as far as both are ^> 1). Therefore, the crossing of Q for 
two different lattice sizes is an evaluation of the critical 
temperature; furthemore, the previous equation can be 
seen as an iteration, in the renormalization group (RG) 
sense. Thus, information on the exponent v can also 
be accessed, as well as other information obtained from 
RG procedures (universality classes, crossover phenom- 
ena, first order phase transitions, etc). For example, for 
the Ising model in zero magnetic field, one apropriate 
function is: 



t = ( sign 




(4) 



where (...) denotes a canonical average, bottom means 
all spins at the bottom plane(linc) of a three(two)- 
dimensional lattice, top stands for all spins at the 
top plane(linc) of a three(two)-dimensional lattice, and 
sign(x) = —1,0 or 1 if x < 0, x — or x > 0, respec- 
tively. For an infinite system, r = 1 for T < T c and 
t = for T > T c . Replacing r for Q in Eq. §, we obtain 
a RSRG equation which connects the parameters K in 
the original lattice and K' in the renormalized (smaller) 
lattice. For a more detailed discussion on the application 
of the FSSRG to the Ising model, see Ref. @. Note, 
however, that the only approximation comes from the 
finite size of the lattices; in fact, even when small lat- 
tices are used, the quantitative results are precise, and 
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this can be understood if one realizes that the FSSRG 
is a generalization, for completely finite clusters, of the 
phenomenological RG developed by Nightingale (l5) . 

In using the FSSRG approach in the study of quantum 
systems, we expect to take advantage of the fact that 
non-commutative aspects of the Hamiltonian are not ap- 
proximated away by the method. Our approach here is 
to use small lattices, which, although not allowing for 
a precise evaluation of the critical parameters (like, for 
instance, those obtained with numerical simulations on 
big lattices), allows for a good qualitative description of 
the critical phenomena involved. The lattices we chose to 
represent the square lattice are depicted in Fig. where 
the right figure is the bigger cluster, with parameters 
K = J/kT and A and the left figure depicts the smaller 
cluster, with parameters K' and A'. The Hamiltonian 
for the original and renormalized clusters are: 

-PH = K[(1- A)(S?Sf + SfS% + SfSf + S v 2 Sl+ 
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and 



- {PHY = K' [(1 - A') (SJSf + S\S V 2 ) 
respectively. 
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FIG. 2. Cells used to implement the renormalization group 
transformation for the anisotropic Heisenberg model on the 
square lattice. 



When applying the FSSRG method to the anisotropic 
Heisenberg model, one has to devise two functions with 
zero anomalous dimension, since it is necessary to renor- 
malize two parameters, namely, A and kT/J. Also, we 
have to distinguish between the antiferromagnetic and 
ferromagnetic models, and two different sets of functions 
were used for each of them. 

For the ferromagnetic model, we chose the following 
quantities: 



t' = (sign(5f)sign(5|)> ; 

rf = (sign(Sf)sign(^)) 
for the smaller cluster and 



(7) 



(8) 



r = (sign^ + S 2 z )sign(S* + S z 4 )) 



V = <sign(Sf + S£)sign(S£ + Sf )) 



(9) 



(10) 



for the bigger cluster. The averages are to be taken with 
respect to the ensemble defined by Eq. || (^) for unprimed 
(primed) quantities. 

For the antiferromagnetic model, Eqs. [f] and ^ remain 
the same, while Eqs. || and |l0| are replaced by: 



<sign(S? 



Sf)sign(- 
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and 



ry = (sign^ - Sf )sign(-Sf + SJ)> , 



(11) 



(12) 



respectively. 

The procedure to calculate these functions is fairly 
straightforward; the final expressions, however, are too 
lengthy and will be omitted. To obtain the required RG 
equations, we impose r' = r and rj' — rj, according to Eq. 
|3|. Approximate values for the critical points are obtained 
from the fixed points of these equations and critical ex- 
ponents are linked to the behavior of the iterations near 
the fixed points. 



III. RESULTS 

A. Finite-size scaling renormalization group 
(FSSRG) 

Our results for the ferromagnetic (F) and antiferro- 
magnetic (AF) models are shown in Fig. ||. The crit- 
ical curve for the F curve is depicted for comparison: 
the accordance with previous results, either approximate 
or exact, is very good. The Ising critical point is lo- 
cated at kT c /\J\ = 2.11, A = 1 (exact values: kT c /\J\ = 
2.269, A = 1); our estimate for the critical exponent v 
for the Ising model is 0.91, while the exact value is 1. 
The universality class for A ^ is the same as for the 
Ising model, which is also consistent with previous re- 
sults. Note that the F and AF Ising models are expected 
to have the same critical exponents and the same modu- 
lus of the critical temperature, and our evaluation agrees 
with these results. 

For the AF model, the phase diagram is qualitatively 
different from its F counterpart: the critical tempera- 
ture reaches zero at a critical value of A, A c , which 
is greater than zero. We find A c = 0.29, which com- 
pares with A c = 0.40 in Ref. g and A c = 0.18 in Rcf. 
||. We also find a reentrant behavior in the critical 
line, which is also present in Rcfs. Nevertheless, 
in Ref. || a second reentrance is observed; the lowest 
temperature we could work with was kT/\J\ = 0.1 and 
we have observed no sign of this second reentrance, nei- 
ther with the FSSRG nor with the bond-moving scheme 
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(to be presented in what follows). In Ref. fl6j, the 
Neel temperature behaves as Tjv ~ l/log(A — A c ) near 
A = A c , while, for the anisotropic ferromagnetic Heisen- 
berg model, T c ~ 1/ log( A) . The supression of long- 
range order at finite temperatures for small values of A 
can be regarded as due to quantum fluctuations. While 
these are not relevant in critical phenomena which take 
place at "high" temperatures, they might be important 
when the critical temperature is low. We expect this to 
be the case for small values of A, and then quantum fluc- 
tuations gain in importance, supressing long-range order. 
Note that we cannot present results at T = 0, since parts 
of our calculation were done numerically and, therefore, 
we are not able to go to very low temperatures. 
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FIG. 3. Critical temperature vs. anisotropy phase-diagram 
for the ferromagnetic (curve F) and anti-ferromagnetic (curve 
AF) anisotropic Heisenberg model. Both lines are atracted to 
the Ising fixed-point, located at kT/\ J\ = 2.11, A = 1. O(D) 
stands for ordered (disordered) phase. The arrow indicates 
the exact critical temperature for the Ising model. 



B. Migdal-Kadanoff approximation 

We have also performed a bond -moving RSRG proce- 
dure to study the AAH model on the square lattice. This 
procedure is equivalent to the one in Ref. |5| applied to 
a different hierarchical lattice; a careful exposition of the 
method is made in Ref. j?| and the bond-moving approx- 
imation is presented in Ref. |16|. 



The phase-diagram for the antiferromagnetic model is 
presented in Fig. ^; note that the qualitative features 
agree with those obtained from other RSRG approaches. 
However, we have not found the second reentrance in the 
critical curve. The Neel temperature Tjv varies as (see 
insert in Fig. 0): 



Tn ~ log (Ac -A)' (13) 

where A c = 0.199 for the Migdal-Kadanoff approxima- 
tion. Again the universality class for the whole critical 
curve is the same as for the Ising model. 




FIG. 4. Neel temperature Tn vs. A phase-diagram for the 
antiferromagnetic anisotropic Heisenberg model on the square 
lattice. The region above (below) the critical line represents 
a(n) desordered (ordered) phase and the insert shows the data 
for small T N , near A c = 0.199. 



We would like to mention that the results from both 
RSRG employed here are in disagreement with spin- wave 
calculations an d numerical simulation |l7],[l8|]. The 
former presents T c = Tn for any value of A and for two 
and three dimensions, while the latter predicts A c = for 
the AAH model. In Ref. |Q , the logarithmic dependence 
of Tn and T c with respect to A — A c is established using 
scaling arguments, but A c = in that paper. Note that 
our calculations cannot be carried out down to T = 0; 
therefore, we cannot study the character of the ground 
state of the AAH model filiT 



IV. SUMMARY 

In summary, we calculate the kT/\J\ vs. A phase- 
diagram for the anisotropic antiferromagnetic Heisenberg 
model on the square lattice. Our results show the pres- 
ence of reentrant behavior on the critical line which sep- 
arates the disordered and ordered phases and a value of 
A > such that the Neel temperature is zero. The en- 
tire critical line is found to belong to the universality 
class of the Ising model. These findings are in agreement 
with previous RSRG procedures (with the exception of 
the second reentrance found in Ref. |J) but not with 
spin-wave calculations and numerical simulation. 

It is clear that the reentrant behavior and a value of 
A c greater than zero are strongly supported by RSRG 
approaches, but the question is not yet settled and more 
work is needed to put these points on firmer grounds. 
Coherent-anomaly methods [ pcfl a nd numerical simula- 
tions with cluster algorithms pl| are possible ways to 
provide a more definite answer to this problem. Work is 
now proceeding along these lines. 
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